function IR = sumar_polar(W,mag)

   [rows,cols] = size(W);
   cy          = round(rows/2);
   cx          = round(cols/2);
     
   radius = min(round(rows/2),round(cols/2))-1;
   angle = 360;
   
   i     = 1;
   suma=0;
   for r=0:radius
      j = 1;
      for a=0:2*pi/angle: (2*pi) - (2*pi/angle)
         if(r>50 && r <90)
           
           x=cy+round(r*sin(a));
           y=cx+round(r*cos(a));
           suma=suma+(W(x,y)*mag(i,j));
              
         end
         
         j = j + 1;
      end
      i = i + 1;
   end
   IR=suma;
   
end